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ABSTRACT 

^ I We present the results from a survey for Extremely Red Objects (EROs) in deep, 

. high resolution optical images taken from the Hubble Space Telescope (HST) Medium 

^ ■ Deep Survey. We have surveyed 35 deep F814W HST/WFPC2 fields in the near- 

0^ I infrared to a typical depth oi K ^ 20. From a total area of 206 arcmin^ and to a 

^ '. limit of is: = 20.0 we identify 224 EROs ((1.14 ± 0.08) arcmin^^) ^-^^^ (^j^^^ - K) ^ 

4.0 and 83 ((0.41 ± 0.05) arcmin^^) with (/gu - K) ^ 5.0. We find that the slope 
of the number counts of the {Isi4 — K) ^ 4.0 EROs flattens beyond K ^ 19, in 
line with results from previous surveys, and the typical colours of the EROs become 
redder beyond the break magnitude. We morphologically classify our ERO sample 
, using visual and quantitative schemes and find that 35% of our sample exhibit clear 

00 ' disk components, 15% are disturbed or irregular, a further 30% are either spheroidal 

\ or compact and the remaining 20% arc unclassifiable. Using a quantitative measure of 

CO . morphology, we find that the ERO morphological distribution evolves across the break 

' in their counts, such that low concentration (disk-like) galaxies decline. We relate the 

fH \ morphological and colour information for our EROs and conclude that those EROs 

Q^' morphologically classified as bulges do indeed possess SEDs consistent with passive 

I ' stellar populations; while EROs with dusty star-forming SEDs are mostly associated 

2 , with disk-like and peculiar galaxies. However, ~ 30% of disk EROs reside in the 

-4—' ' passive region of I/J/K colour-colour space. These could be either genuinely passive 

^ I systems, lower redshift contaminants to the high-z ERO population, or systems with 

• • , composite star- forming and passive SEDs. We use photometric redshifts for our high 

^ ^ ' S /N multicolour photometry and derive redshift distributions in good agreement with 

\ spectroscopic work of somewhat brighter ERO examples. 

^ . Key words: cosmology: observations, galaxies: bulges, galaxies: peculiar, galaxies: 

■ - - ' starburst, galaxies: evolution, infrared: galaxies 
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1 INTRODUCTION 

Broadband photometric selection is an increasi ngly common 
technique to identify high redshift galax ies iSteidel et alJ 
ll999l:lTotani et all20oil : lFranx et al.l2003l) . One particularly 
powerful combination is to use optical and near-infrared 
passbands, typically Rov I and K, to identify galaxies with a 
strong decline in their spectral energy distributions (SED) at 
wavelengths around l^m. These sources display extremely 
red optical-near-infrared colours, e.g. (7 — K) ^ 4.0, and 
as such are termed Extremely Red Objects (EROs). This 
strong spectral decline could either result from a break in 
the SED or the suppression of bluer light due to strong dust 



* Based on observations from the NASA/ESA Hubble Space Tele- 
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absorption. As we look at more distant galaxies, the first 
spectral break strong enough to produce such a red colour 
is the 4000A break, which falls between the / and K pass- 
bands at z; > 1. Progressively bluer features, such as the 
Balmer break (3625 A) and Lyman-limit (912 A), will fall 
in the relevant wavelength range for higher redshift sources, 
although these will typically be fainter than those systems 
selected on the basis of the 4000A break. Similarly, the dusty 
systems within the ERO population are also expected to lie 
at high redshifts, due to the strong wavelength dependence 
of dust reddening combined with the redshifting of UV light 
into the optical passbands. ERO surveys thus enable us to 
identify both the most and least active galaxies at 2: ^ 1, 
and hence pro vide a powerful pro be of obscured star forma- 
tion at z ^ 1 iSmail et al.ll2002bl hereafter S02) as well as 
the nature of the evolved descendants of some of the earli- 
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est phases of sta r formation in the Universe jPunlop et alJ 
ll996l:IZepj[l997lV 

Unfortunately, the mixed nature of the ERO popu- 
lation, compounded by the s trong angular clu stering of 
one or both sub-populations IPad di et al.ll200(l l. has led 
to a vigorous debate about the relative importance of the 
evolv ed and active comp onents in the overall ERO popula- 
tion ijCimatti et alJl2002t S02). 

The deepest ERO surveys undertaken show a rapidly 
increasing surface density of sources as a function of mag- 
nitude limit down to -ft' ~ 19. Fainter than this, the in- 
crease in the numb ers of EROs with apparent magnitud e 
slows dramatically (|Smith et alJ l2002al : iFirth et al.ll2002l) . 
This transition is very abrupt and suggests that there may 
be a profound change in the nature of the EROs fainter than 
K ~ 19, perhaps associated with a variation in the relative 
propo rtion of passive and active systems at fainter magni- 
tudes iSmith et al.ll20o'2al : S02) . To investigate this possibil- 
ity we need an observational test to distinguish between the 
different sub-classes of EROs. There are three approaches 
which have been used to attempt to differentiate passive 
from dusty/active EROs: photometric tests, dust-insensitive 
star formation indicators and morphological classification. 

The photometric classificat ion of EROs was popularised 
by IPozzetti fc Mannuccil (1200011 , using J-band photometry 
to separate the two sub-classes on the (/ — K)-(J — K) 
plane. Subsequent application of this test to various sam- 
ples has led to the conclusion that the population at /sT ^ 20 
is roughly equal ly split between the dusty /active and pas- 
sive sub-classes l|Manmicci et al.ll2 002l: IVa iizella et 80120021 : 
ISmith et all l2002al : see also ICimatti et al.ll2002ll . Unfortu- 
nately, without an independent test of this classification 
scheme it is hard to know how reliable these conclusions 
are. 

Searches for millimetric and radio emission from sam- 
ples of EROs have also been used to attempt to estimate 
the proportion of dust-obscured star-forming EROs in the 
overall population llMohan et al.ll2002l: IWehner et "811120021: 
S02). However, these star- formation indicators, while being 
insensitive to dust obscuration in the galaxies, are only ca- 
pable of detecting the most vigorously star-forming systems 
at 2 ^ 1. As a resul t they provide only a lower-limit of 
> 10% l|Mohan et alj E'002:l and > 30% (S02) on the pro- 
portion of star-forming systems in the ERO population. S02 
also compare the colours of their radio-selected star-forming 
EROs with those expected for dusty star-forming and pas- 
sive galaxies and find that most fall on or within the model 
classification b oundary for dusty galaxies, providing some 
support for the lPozzetti fc Mannucci i2000l) scheme. 

The final observable which has been used to differentiate 
between evolved and active ERO populations is their mor- 
phologies. It is commonly assumed that t he passive EROs 
will have elliptical or at least early-type iMoriondo et alj 
I2OOOI: ISmith et al.ll2002bll morphologies, while star-forming 
galaxies are expected to exhibit eithe r a disk-like, dis turbed 
or obviously interacting morphology llDe^eTailiggi). Mor- 
phological studies of very red galaxies h ave indeed shown 
that some have elliptical m orphologies jMenanteau et al] 
I1999I: iTreu fc StiaveUil ll99^. However, results on statisti- 
cal sample s of EROs have produ ced more mixed c onclu - 
sions, with iMoriondo et al.l l^2000^ and iRoche et all 
claiming a high early-type fraction and more recent work 



bv lSmith et ai] and lYan fc ThompsonI i2003l here- 

after YT03) arguing for a more equal mix of early-type and 
later-type/disturbed systems. By comparing this classifica- 
tion scheme with those described above it is possible to test 
the reliability of the diffe rent approaches for a large sample 
of EROs (c.f. lSmith e~l. 2002a). 

We highlight the study of YT03 which has very similar 
goals to those described in this paper. Their work utilised 
high resolution FST WFPC2 F814W-band imaging to mor- 
phologically classify a sample of 115 (/814 — K) ^ 4.0 EROs 
selected from ground-based A's-band imaging to a 5-(t depth 
of Ks S 18.7. The authors visually classified their objects 
and found a mix of around 30:64:6 for bulge dominated, 
disk dominated and unclassified classes, with around 17% of 
the sample showing signs of merger/strong interaction. They 
showed their visual estimates to be in generally good agree- 
ment with an automated bulge-|-disk decomposition tech- 
nique. The key improvement of our study over that of YT03 
is that our near-infrared imaging is sufficiently deep to probe 
beyond the break in the ERO counts at ~ 19. We also im- 
prove upon their work through the addition of multicolour 
photometry to test the different techniques for subclassifi- 
cation of the ERO population. 

We discuss our near-infrared observations of a selected 
sample of deep Hubble Space Telescope WFPC2 images in 
§2, present our analysis of these in §3 and the results which 
this provides are presented and discussed in §4, before giving 
our conclusions in §5. Throughout we assume a cosmology 
with Q.m = 0.3, f2A = 0.7 and Ho = 70kms"^ Mpc"\ 



2 OBSERVATIONS AND DATA REDUCTION 

A common misconception is that the key to a successful 
ERO survey is to obtain deep NIR imaging - in fact, the 
observationally most demanding aspect is achieving the nec- 
essary depth in the optical to identify that a galaxy is an 
ERO. For this reason, we have chosen to concentrate our 
survey on fields for which deep, high-quality optical imag- 
ing alrea dy exists. These field s come from the Medium Deep 
Survey fCriflit hs etaP Il994l) which consists of over 500 
deep HST /W¥¥C2 images of intermediate/high-Galactic 
latitude blank fields. We selected a subsample of ~ 100 fields 
from the parallel and pointed samples from the survey. These 
were selected by requiring only that the total exposure time 
in the F814W filter be greater than 4.0ks, that the primary 
target of any parallel observations was not a galaxy cluster 
or that, if the observations were pointed, that they were not 
targeting an extragalactic source (to ensure representative 
extragalactic regions) and that they be suitably placed for 
northern hemisphere follow up (Dec. > —10°). These then 
represent high-resolution (0.1 arcsec FWHM) and very deep 
{him ~ 25-26) images of random areas of the extragalactic 
sky. 



2.1 Near infrared observations 

Near infrared imaging was obtained u sing the Isaac Newton 
Group Red Imaging Device CINGRID. IPackham et 'al]l2003h 
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on the 4.2m William Herschel Telescope (WHT)\ INGRID 
IS a 1024^ HAWAII-2 array at the bent Cassegrain focus of 
the WHT giving a 0.238 arcsec pixel"^ scale and a 4.1x4.1 
arcmin field. 

Observations were made over 12 nights, 2000 November 
13-15, 2000 December 9-11 and 2001 May 1-6. One night 
was lost due to the instrument being unavailable. A total of 
55 fields were observed in the ifs-band and 30 of these in 
the J-band. Four of the nights were non-photometric, and 
we are in the process of obtaining calibration data for these. 
For the rest of this paper, we consider the bulk of the Kb- 
band data, 35 fields covering 206 arcmin^ [HST + _ft'-band 
coincident imaging) which are well calibrated. A summary 
of these fields is given in Tabled 

Each MDS field was imaged using a 9 point dither pat- 
tern with exposures of around 60s in Ks and J exposures of 
around 120s. Each image was made up of 4 x 15s expo- 
sures coadded in hardware. The 9 point dither pattern itself 
was moved around the target field, in order to avoid bright 
objects falling on the same pixels. Total integration times 
are typically 2.8ks and are listed in Tabled 

The INGRID data were reduced us- 
ing a custom written pipeline, available from 
|http : //star-www . dur . ac.uk/~ dgg/ipipe/ The pipeline 
uses standard iraf'^ and starlink routines, and the 
reduction procedure is as follows. 

2.1.1 Initial processing 

The first step in the processing was to construct a bad pixel 
mask and apply a correction for an offset in the exposure 
time present in early versions of the software controller^. 
Dark subtraction was found to be unnecessary for most 
frames, but a master dark frame was subtracted if resid- 
ual structure was seen. The data were flatfielded and sky 
subtracted using the in-field c hopping, or moving flatfield 
technique of lCowie et alj lll990fl . using a ruiming median of 
8 temporally adjacent frames. 

2.1.2 Mosaicking 

Image registration and mosaicking was performed entirely 
using the starlink software CCDPACK. The first image of 
each target was taken as the reference and the relative offsets 
of the others calculated. This was achieved using the tasks 
FINDOBJ, FINDOFF and REGISTER to find objects in common 
between frames and calculate the offsets to sub-pixel accu- 
racy. The images were then geometrically transformed with 
subpixel shifts and bilinear interpolation to conserve flux 
using TRANNDF, and finally combined using a 3-a clipped 
mean within the makemos task. 

^ Based on observations made with the William Herschel Tele- 
scope operated on the island of La Palma by the Isaac Newton 
Group in the Spanish Observatorio del Roque de los Muchachos 
of the Institute de Astrofisica Canarias. 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatory which is operated by AURA Inc. under contract with the 

N SF. 

^ |http : / /www . ing . iac . es/Astronomy/instrument s/lngrid/ 1 
ingrid_timing.html 
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The mosaic so produced is regarded as a first pass mo- 
saic. In order to improve sky subtraction, a mask is made 
for all the astronomical objects in this image. SExtrac- 
TOR v2.2.2 jBertin fc Arnoutj|l996li was used to detect and 
mask objects by making a CHECKIMAGE with the OB- 
JECTS option. 

The sky subtraction step is then repeated by first ap- 
plying the de-registered mask to the images going into the 
local flatfield. In this way, objects too faint to be found in 
the individual exposures which would otherwise bias the es- 
timate of the sky level and lead to oversubtraction of the sky 
can be successfully rejected (M. Currie, priv. comm.). The 
mosaicking step is then repeated with these second pass sky 
subtracted images to make the final mosaic. 

In addition to the image mosaic, an exposure map is 
constructed by summing the number of images going into 
each pixel. This is used to deal with the uneven noise prop- 
erties of the mosaics. An image with uniform pixel to pixel 
noise is generated by multiplying the mosaic by the square 
root of the exposure map. However, for the current dataset, 
the field of view of INGRID is very well matched to the 
WFPC2, allowing the latter to be entirely contained within 
the area of the former for any roll angle. Thus, no region of 
the WFPC2 field sees less than 100% of the near-infrared 
(NIR) exposure time, except the three pairs of fields suf- 
fixed 'XA' and 'XB'. These are fields where two overlapping 
MDS pointings exist, and the INGRID pointing lies midway 
between. Only NIR data with more than 50% of the total 
exposure time is considered in these three fields. 



2.1.3 Non-photometric data correction 

In the discussion of the mosaicking process above, no correc- 
tion has been made to the photometric zeropoints for each 
individual frame going into the mosaic. This is because the 
extinction in the NIR is low {<, 0.08 magnitudes/airmass) 
and the fields were typically observed for around 2.8ks at 
low airmass, therefore the change due to the extinction vari- 
ation from image to image is negligible. This is not the case 
for non-photometric data. However, it is straightforward to 
correct all the exposures of one field to the same zeropoint, 
as the interval between exposures is short, and the dither 
pattern is small, so there are many objects in common be- 
tween frames. 

To correct for non-photometric data, the first frame of a 
field was again taken as the reference image. The first pass 
mosaic was made as before, but when SExtractor was 
run on it to make the object mask, a catalogue of the 20 
brightest objects in the field was output. After the second 
pass sky subtraction, aperture photometry was performed 
on these 20 objects and the optimal photometric zeropoint 
shift between each exposure and the reference calculated 
using the linfit task in iraf. Typical shifts were ~ 0.1 mag. 
The validity of this approach is verified through comparison 
with repeat observations made in photometric conditions, 
and with external photometry, described next. 

2.1.4 Photometric and astrometric calibration 

Calibration on to the UKIRT Tf-band syste m was performed 
using UKIRT Faint Standard (FS) stars jHawarden et alJ 
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Table 1. Table of observations. " — Field names are derived from MDS catalogue names. Fields suffixed XA and XB represent two 
overlapping HST exposures, with one corresponding /S'-band exposure targeted at the midpoint of the two; ^ - Target name is target of 
HST observation for pointed fields; PAR indicates a pure parallel field. 
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l20nifl . Several such stars were observed at the start and 
end of each night. The accuracy of the calibration was 
checked against the 2MASS All-Sky Point Source Cata- 
logue*. Around 60% (including non-photometric nights) of 
the fields possessed sufficient bright objects in common with 
2MASS to obtain an accurate test of the calibration. The 
minimum number of objects used was five, and more typ- 
ically 10 - 20. Where the data overlapped, the agreement 
was found to be typically better than 0.05 mags. This also 
showed that zero point on non-photometric nights did not 



differ from the nominal zero point estimated from the stan- 
dard stars by more than ~ 0.3 mags. We stress that we do 
not use data where the absolute calibration is more uncer- 
tain than ~0.1 mags. All passba nds are corrected for G alac- 
tic reddening using the maps of ISchleeel et all lll99dl . As- 
trometry was performed against the USNO-A2.0 astro metric 
catal ogue l)Monetll998h using the package wcstools (|Minkl 
Il997tl to automatically calculate the World Coordinate Sys- 
tem (WCS) of the image, and found to be accurate to better 
than 1 arcsec. 



This publication makes use of data products from the Two Mi- 
cron All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis Cen- 
ter/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 



2.2 HST imaging 

The HST images were retrieved from the ST-ECF Associa- 
tion archive. The archive takes all available associated data 
defined by the pointing of the telescope, performs an on-the- 
fly recalibration (OTFR) and co-adds these chip by chip. 
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The only additional processing necessary was to combine 
the four WFPC2 chips into a single mosaic. In so doing, the 
smaller pixel scale of the PC chip is resampled to that of the 
three WF chips. This was done using the iraf stsdas task 
WMOSAIC, which also corrects for the geometric distortion 
within the instrument. Photometric zeropoints were taken 
from the headers produced by the OTFR and are cal i brate d 
to the Vega system as described in lHoltzman et al] il995l) . 
The F814W Vega system is very close to the Cousins /-band 
and shall be referred to as i814 henceforth. 



2.3 Image alignment and photometric catalogue 

The seeing in the ground-based data was automatically de- 
termined by profile fitting to stellar objects in each frame. 
For each field, the HST data were then matched using Gaus- 
sian convolution to the seeing of the passband with the poor- 
est seeing. All passbands were aligned to the _ft's-band image 
and resampled to the INGRID pixel scale^, again using the 
STARLINK package CCDPACK to perform automated registra- 
tion. The field of INGRID is astrometrically very fiat, as 
is the WFPC2 field after geometric correction with WMO- 
SAIC; thus the residual offsets between objects in the dif- 
ferent passbands were typically at the level of one pixel. 
Although the internal astrometry of the WFPC2 device is 
very good, the external accuracy is somewhat lower, thus the 
WCS that we adopt in the catalogues is that of the Ks-hanA 
image. 

Photometry is performed on coaligned images after con- 
volution to the same PSF. These shall be referred to as the 
convolved images. Also, in order to exploit the exquisite reso- 
lution of the HST, we will wish to use the images prior to this 
transformation. We shall refer to these images as the non- 
resampled images. The CCDPACK software allows the map- 
ping between our different coordinate systems to be readily 
stored, and a provides convenient way to transform between 
them. 

SExTRACTOR was used to detect objects brighter than 
the fiK ~ 21.2 mag. arcsec"^ isophote, with at least 4 con- 
nected pixels (0.23 arcsec^), after filtering with a Gaussian 
with a FWHM of 4 pixels (0.95 arcsec). This limit is de- 
termined from our shallowest fields, in order to ensure high 
uniformity in the object catalogues. Object detection was 
run on the uniform noise mosaic, described above. Photom- 
etry is performed by running SExtractor in double im- 
age mode, using the detection image to define the location 
of photometric apertures and measuring magnitudes from 
each convolved passband in turn (note that if-band mag- 
nitudes were measured from the A's-band image, not from 
the iiTs-band uniform noise image) . An aperture diameter of 
2.5 times the FWHM of the seeing was used. This ensure s 
high signal to noise within the aperture ijLillv et al.lll99ll) . 
and is small enough (typically around 2 arcsec) to avoid 
contamination from neighbouring objects, due to the good 
seeing of our data. We verified that our choice of photomet- 
ric aperture is robust, by also examining colours defined by 
a fixed 2 arcsec aperture, a 4 arcsec aperture, and 'total' 
colours measured using the SExtractor BESTJVIAG in 

^ Hereafter, pixel shall refer to this common pixel scale of 0.238 
arcsec unless explicitly stated otherwise. 



each band (where the BEST_MAG aperture is defined from 
the ifgimage). It was found that all these measures are con- 
sistent, but that our adopted approach typically yields the 
colour with the smallest error. 

Total magnitudes are measured for the K-hand us- 
ing the SExtractor BEST_MAG. 5-a limiting magnitudes 
within these apertures were calculated using SExtractor's 
estimate of the signal to noise for each object detected (Ta- 
ble 4) . The adopted 5-cr limiting aperture magnitude for the 
survey is if ^ 20.0, and only objects brighter than this are 
considered. 

To assess completeness, a set of simulations was run first 
by inserting artificial point sources and then artificial ellip- 
tical and disk galaxies (of the typical size of faint galaxies 
in the survey) into our Ks-hand images. Our detection pro- 
cedure was then repeated on these. The typical complete- 
ness for point sources is greater than 80% ai K — 20.0, 
and a similar completeness level is reached around 0.3 mag 
brighter for objects with de Vaucouleurs profiles. The sim- 
ulations show we are 80% complete for compact sources at 
K=20.0; while they also show that the sample is 100% com- 
plete for point sources around K=19.4, de Vaucouleurs pro- 
files at K=19.1 and the intermediate case of exponential disk 
profiles is 100% complete around K=19.2. Since our primary 
aim is to interpret the morphologies of EROs in the context 
of our survey and previous imaging surveys, which are typi- 
cally selected to a limiting 5-a depth, we apply this selection 
to our sample. A discussion of the completeness for differ- 
ent morphological types will be considered in future work, 
when comparing our data with galaxy formation models. In 
order to assess the spurious detection rate, the detection 
images were multiplied by —1 and the detection and mea- 
surement procedure repeated. No false detections brighter 
than K = 20.0 were found. 



3 ANALYSIS 

Our analysis will proceed by using our K-selected galaxy 
survey to isolate Extremely Red Objects based on their 
(^814 — K) colours. We will use the superlative resolu- 
tion of the HST imaging to morphologically classify these 
EROs. Both qualitative visual classification and automated 
machine-based classification schemes will be used. We shall 
then compare our photometry, number counts and morpho- 
logical fractions with previous work. 

3.1 ERO sample selection 

EROs were identified by selecting all objects match- 
ing the common definition (e.g. S02, IWehner et alj |2002| . 
iRoche et alj2003l. YT03) of an ERO of {Isi4-K) ^ 4.0. Our 
goal is to determine the morphological mix of such EROs to 
K~20. The astrophysical motivation behind this colour se- 
lection is described in detail in Pozzctti & Mammcci ( 20(^ 
and references therein. Briefiy, this colour cut corresponds 
to the expected colours for a z~l passively-evolving, old 
stellar population, but is much redder than the expected 
colours of normal field galaxies; only high-redshift, heavily 
dust-reddened, massive starbursts can have similar colours. 
Therefore, this colour cut isolates the most and least actively 
star- forming galaxies at z^l. 
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Fig. shows the colour-magnitude diagram for the 
whole -ft'-selected survey. Our ERO catalogue was visually 
inspected and a small number of sources were rejected as 
being false detections, either being saturated stars or ob- 
jects lying in the periphery of bright stars. No other false 
detections were found. We emphasize that for our entire sur- 
vey, the NIR data used is highly uniform across the entire 
WFPC2 field of view and does not suffer from reduced sig- 
nal to noise at the edges, as can be the case for surveys 
undertaken with smaller NIR imagers (e.g. YT03) . We class 
objects with a SExtractor CLASS_STAR ^ 0.97 as galax- 
ies. Only one of our EROs would be classed as stellar by this 
definition (ERO082). We include this object for complete- 
ness (it may still be extragalactic, e.g. AGN, and will be 
morphologically classified from WFPC2 imaging later) but 
note that its inclusion has a negligible effect on our results. 

We identify 224 sources to 20.0 with {I - K) ^ 4.0 
and 83 with (/ — K) ^ 5.0 across the 206 arcmin^ of our 
survey. 

3.2 Morphological analysis 

Using the exquisite resolution of the WFPC2 we can mor- 
phologically classify the EROs in our survey. This classifica- 
tion will proceed in two ways: we will use a scheme based on 
visual inspection of the images and also a quantitative sys- 
tem using machine-based measurements of the central con- 
centration of each object. Each method has its own advan- 
tages and disadvantages: visual classification is somewhat 
subjective, but is better suited to deal with unusual objects, 
not well-fit by a model; whereas automated classification 
should be reproducible for a given set of input parameters, 
but unusual objects may lead to catastrophic misclassifica- 
tions - e.g. for a merging system, the software might only fit 
one subcomponent, considering it to be an isolated system. 

For the morphological analyses, the non-resampled 
F814W images were taken, and residual cosmic rays were re- 
jected using the iraf task CRAVERAGE. Next the images were 
rebinned 2x2 to increase the per pixel surface brightness 
sensitivity. The reduction in angular resolution from this 
procedure does not adversely affect the morphological classi- 
fications, but aids in classifying the lowest surface brightness 
systems, which would otherwise be unclassifiable. 

3.2.1 Visual classification 

All the EROs were visually classified by one of us (IRS) 
following the scheme in Table |5| This scheme was devised to 
isolate the broad classes of objects identified on a first-pass 
through the dataset. 

This procedure used the same display software em- 
ployed by the MORP HS group when cla ssifying faint galax- 
ies in distant clusters fSm ail et al]|l99i1l . The software dis- 
plays two images at difl^erent stretches of a 10" x 10" region 
around the ERO from the non-resampled WFPC2 images, 
after they are smoothed with a 0.2" FWHM Gaussian to 
reduce the shot noise. This was used in conjunction with a 
hard-copy of each thumbnail image displaying the Ks-haxiA 
image of the ERO contoured over the F814W image to allow 
isolated i814 -band components within a single ERO to be 
easily identified (e.g. Fig. II 111 . 



Table 2. Visual classification scheme for the EROs, based on the 
WFPC2 F814W imaging. Classes are motivated by the properties 
of the sample and we bin these to describe disk-like (2, 5, 6) and 
bulge-like (3 and 7) systems. 



Class 


Description 


# 


% 





Blank or too faint to classify 


43 


19.20 


1 


Compact (small and peaked) 


32 


14.29 


2 


Compact, disk 


36 


16.07 


3 


Compact, symmetrical 


26 


11.61 


4 


LSB, disturbed 


1 


0.45 


5 


LSB 


6 


2.68 


6 


Obvious disk 


42 


18.75 


7 


Spheroidal 


7 


3.12 


8 


Merger 


13 


5.80 


9 


Amorphous 


18 


8.04 



We give the proportions of difi'erent morphological 
classes in our full K^20 sample in Table |5| In a testament 
to the depth of the HST imaging, less than 20% of these 
very red and very faint galaxies were unclassifiable in the 
WFPC2 images. For the remaining sources, 40% are com- 
pact (although these can still exhibit weak morphological 
features), a further 20% are relatively well- resolved, disk 
galaxies and the final 20% comprise a mix of low surface 
brightness (LSB), mergers, amorphous or clearly spheroidal 
galaxies. As we have stressed, these morphological classes 
were motivated by the properties of the ERO sample and 
hence we will need to combine and convert them to trans- 
late them into astrophysically interesting classes. 

3.2.2 Automated classification 

Given the modest S/N of typical EROs in our WFPC2 
imaging, we have chosen in our quantitative morphological 
analysis to concentrate on two simple measures of the light 
distribution within the EROs. These are t he concentration 
index, C, an d mean surface brightness, /i iAbraham et alJ 
Il994ll996allbl) . Measurements of C for the ERO sample were 
made using a modified version of SExtractor to calcu- 
late C — Fo.a/Fi.o where is the integrated flux within 
an elliptical aperture which contains a fraction x of the 
total isophotal area, and x=0.3 is the standard definition 
(Abraham ct al. 199^. The error on this measurement was 
estimated from simulations. Grids of artificial galaxy images 
were constructed using the iraf ARTDATA package and in- 
serted into random blank sky regions of the non-resampled 
F814W image, and the scatter in the measurement of C for 
the same input profiles recorded as a function of apparent 
brightness of the galaxy. The measured variation in C with 
apparent brightness of the galaxy (for a fixed light profile) 
in these simulations also allows us to investigate the evolu- 
tion in the morphological mix of the EROs as a function of 
magnitude in a quantitative manner. 

3.3 External comparisons 

3.3.1 K-band galaxy counts 

As a check of our AT-band photometry, in Fig. |5| we com- 
pare our galaxy number counts (for all K-selected galaxies. 
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Figure 1. Colour-Magnitude diagrams for all galaxies in the survey (left panel) and for EROs (right panel) defined by (/814 — K) ^ 4.0. 
Error bars are omitted from the left plot for clarity and the broken line illustrates the (/814 — K) ^ 4.0 ERO selection criterion. EROs 
are highlighted with bold points. Note the rapid increase in the proportion of the galaxy population in the ERO class at > 18 (4%), 
followed by a more constant distribution fainter than K ^ 19 where the EROs contribute 12% of the total population. 
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Figure 2. Differential X-band galaxy counts in half-magnitude 
intervals, corrected for incompleteness. Error bars simply assume 
Poisson statistics and are thus likely to be underestimates of the 
true errors. Selected number counts from the literature are also 
plotted. Our counts are in good agreement with other surveys. 

irrespective of colour) with a selection from the literature. 
Our counts are in good agreement, showing a steep rise with 
gradient a = 0.34 ± 0.06 to K'~20, with an incompleteness- 
corrected surface density of (7.5 ± 0.2) arcmin"^ at K=20. 

3.3.2 ERO number counts 

Firstly we consider the number counts of (Isu — K) ^ 4.0 se- 
lected EROs in our survey. We compare our number counts 
to other datasets which use the same (/8i4 — K) selection to 
remove any concerns about differences in the samples due 
to the photometric selection (although see t|4.71 . This com- 
parison is shown in the left panel of Fig. |3 Our cumulative 
ERO counts appear to show a break around K ^ 19 and 
so we fit them using two power laws (log^Q N{Ki < K < 
K2) oc aK), giving slopes of: a = 0.88 ± 0.09 for K < 19.0 
and a = 0.42 ± 0.19 for 19.0 < ^ 20.0. Also shown in 



the figure are the counts of (/814 — K) ^ 4.0 EROs from S02 
(priv. comm.), Wchncr ct al. (2002), Roche ct al. (2003) and 
YT03. 

We se e that our ERO c ounts are in excellent agree- 
ment with lRoche et all (|200^, based on publicly available 
ESO imaging of the CDFS/GOODS. The S02 counts are 
marginally lower than the other data, but note that the er- 
ror bars are purely Poissonian errors an d likely to underes - 
timate the true field-to-field variation JPaddi et alJl2000l) . 
At the faint end t heir number density agrees well with both 
ours and those of iRoche et"ai] (12003) . The counts of YT03 
appear marginally higher but not significantly so. 

We find an incompleteness-corrected surface density for 
(/814 - K) 4.0 objects of (1.14 ± 0.08) arcmin"^ to K ^ 
20.0, representing 17% of the total galaxy population to this 
depth. This is in reasonable agreement with the value of 
(0.94 ± 0.11) arcmin"^ from S02 (using th e same selection 
limits as ours, priv. comm. ). IWehner et al] (1200211 also find 
comparable numbers, but with large field to field scatter 
between their three fields, each of which uses a cluster lens 
to increase the sensitivity of their observations. The data 
shown in the plot excludes one of their three cluster fields 
(A2390) which seems particularly overdense in EROs. This 
further emphasizes the impact of field to field variations. The 
surface density to the magnitude limit of the YT03 survey 
also appears in good agreement with our estimate, and is 
consistent with the S02 point. 

Next we consider a redder subsample selected by (/814 — 
K) ^ 5.0. We find a surface density of such objects of 
(0.41 ± 0.05) arcmin"^ or 5% of the population to K ^ 20. 
The number counts of (/814 — K) ^ 5.0 objects are ad- 
equately described by a single power law with a slope of 
a = 0.70 ± 0.02. Here we can only compare with the 
(I — K) ^ 5.0 counts of S02 (priv. comm.) who find a surface 
density of (0.18±0.05) arcmin~^. Again we note that the er- 
ror estimates do not include a contribution from clustering. 
We also note that this extreme colour cut is on the rapidly 
falling tail of the colour distribution (encompassing ~5% of 
the population). Slight photometric offsets between differ- 
ent surveys leading to different colour measurements could 
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cause a large number of objects to be included or excluded. 
For example, if our two surveys measure colours systemat- 
ically different by 0.2 mags (a reasonable calibration offset 
between ground based surveys), this is enough to allow our 
Poissonian error bars to overlap. 

3.3.3 Comparison with YT03 photometry 

We compare our photometric measurements with those of 
YT03. We have f 5 MDS fields in common, from our sample 
of 35 and their sample of 77 MDS fields. Firstly we take 
the FRO catalogue of YT03 (their Table 2) and search for 
objects which lie in our survey fields (Table 0. We find 28 
such objects (only 17 of which do we classify as EROs). Next 
we search our galaxy catalogues for entries within l"of the 
coordinates given by YT03. We match 26/28 objects. The 
remaining two (their u2vl6#108 and u2vl9#30) lie between 
the WFPC2 chips in our images and were thus masked from 
our analysis. The ST-ECF associations archive only includes 
data taken at very similar roll angles, so it is possible that 
YT03 included additional imaging not used in our analy- 
sis. Comparison of our F814W exposure times for the fields 
in question (U2V16 and U2V19) show identical exposure 
lengths to theirs. Therefore we believe that their inclusion 
of these objects in their catalogue is spurious. 

We compare our K-hand total magnitudes with their 
large-aperture ('total') magnitudes for these 26 galaxies and 
plot the results in Fig. |I| For 10/26 (38%) our estimates 
agree to within 1-cr, however, for 13/26 (50%) EROs YT03 
measure considerably brighter K-hand magnitudes, by 0.5- 
1 mag. This offset thus results in YT03 measuring redder 
(-^814 — K) colours for these galaxies and hence classifying 
them as EROs, whereas in our data they are bluer than 
(IgiA — K) ^ 4.0. We visually examined the extreme outliers 
and can find no obvious reason for the discrepancy. For ex- 
ample, one such outlier with AK = 1.21 lies on the same 
WFPC2 chip as two other EROs which we measure to be 
the same brightness as YT03 to within 0.4 magnitudes. For 
these 26 galaxies, we measure a median offset of 0.26 ±0.07. 
This slight photometric offset would bring thei r ERO num- 
ber co unts into closer agreement with ours and lRoche et alJ 
l)2003() . We emphasize that this comparison only uses data 
taken under photometric conditions and so cannot be af- 
fected by corrections for non-photometric images. We at- 
tribute the disagreement to the quality of the YT03 data, 
which is known to be shallower, more inhomogeneous and 
taken in poorer seeing conditions. For the case illustrated 
above, the seeing FWHM in their data is 1.7", compared 
with ours which is around 1.0". The poorer seeing means 
that a larger photometric aperture must be used, leading to 
a noisier measurement of the magnitude. 

3.3.4 Comparison with YT03 visual morphologies 

We directly compare our visual classifications with the clas- 
sifications of YT03, for the 17 EROs in common between our 
samples. Though our classification schemes are different, we 
can bin both samples coarsely into disks, bulges and other. 
Our classes 2, 5 and 6 are disk dominated classes and broadly 
equivalent to their D, D-l-B and ID classes. We find 11 ob- 
jects in common in these three classes. For spheroidal des- 
ignations, we compare our classes 3 and 7 with their B and 
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Figure 4. The comparison between our i^-band photometry and 
that of YT03 for galaxies in common. We plot the difference be- 
tween their large aperture magnitude and our total magnitude. 
Even though our total magnitude should encompass more light 
than their aperture magnitude, there is evidence that their K- 
band magnitudes are systematically brighter, by around 0.3 mag 
(indicated by broken line). See text for discussion. 

BD. Here we find that one object we classify as spheroidal, 
YT03 class as a disk. Furthermore, two EROs which they 
class as bulges, we classify as disks. From our merger cate- 
gory (which has no equivalent in their work), we have two 
objects, one of which YT03 class as a bulge, the other as 
a disk; and we both agree that a final object is unidenti- 
fiable. Hence, through a very broad comparison with the 
YT03 classifications, we agree on 12 classifications, disagree 
on 3 and a further 2 (our merger category) we are unable to 
compare. 

Next we examine the distributions of morphologies 
within our ERO samples. We restrict our catalogue to K ^ 
18.7, the median 5-a limiting magnitude of YT03, in order 
to compare the morphological mixes between the two sur- 
veys. From their visual classifications, YT03 found 35 ± 5% 
of their EROs to be bulge-like and 64 ± 7% to be disk-like 
with 6% unclassifiable (where the errors are simply Poisso- 
nian). Our visual classifications divided coarsely into bulge- 
like or disk- like classes (Table|2l yield 24 ±6% bulge-like and 
54± 10% disk-like with 8% unclassifiable. Note that here the 
14% (7/50) of our classifications which do not fit into the 
simple bulge-like or disk-like (e.g. merger, amorphous) are 
not included. Hence, both our surveys find approximately 
twice as many disks as bulges at these magnitudes. 

3.3.5 Comparison with MDS automated morphologies 

A more detailed quantitative an alysis of galaxy morphol- 
ogy in the MDS was performed bvlRatnatunga et alJ (Il994l . 
1999). In outline, they applied a maximum likelihood fitting 
technique to the surface brightness profiles of MDS galaxies. 
The models they fitted were axisymmetric de Vaucouleurs 
and exponential disk profiles. A galaxy best fit by the for- 
mer profile was classed as a bulge and the latter profile as a 
disk. If some combination of the two gave a better fit then 
the relative contribution of the two components was out- 
put in the form of a bulge to total luminosity (Ls/hTotai). 
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Figure 3. Surface density of EROs. Left panel shows cumulative number counts for the full ERO sample selected with (/gi4 — K) ^ 4.0. 
Overplotted are other ERO counts from the literature, with their selection criteria listed. Right panel shows the redder subsample selected 
with (/gi4 — K) ^ 5.0 compared with the redder subsample of S02 (priv. comm.) using the same selection criterion. Solid lines are our 
power law fits. 



The completeness limit for classification is drawn at H = 1.8 
(where H is their signal to noise index computed from the 
signal to noise of each pixel in an object) which corresponds 
to / ~ 24.5 for the shallowest fields used here. A H ^ 2.0 (0.5 
mags brighter) is required to class a galaxy as either disk-like 
or bulge-like (otherwise sources are just divided into either 
point or galaxy), and a H ^2.4 is required to fit a bulge-|-disk 
model. 

The reliability of this maximum likelihood estimate 
(MLE) was tested with the sophisticated and now widely- 
used two-dimensional surf ace brightne ss fitting algorithm of 
ISimard (,I99&. GIM2D) bv lSimard et al. (2002). They found 
that the agreement of structural properties such as half-light 
radius and bulge to total light ratios between the classifiers 
was very good down to the limit studied of /814 55 22. A 
small fraction of objects were found to have very different 
bulge fractions as measured by the different classifiers, and 
most of these turned out to be peculiar or interacting sys- 
tems. This emphasizes that automated classifiers are unable 
to deal robustly with unusual galaxies. 

Morphological fractions for our ERO sample based on 
the MDS MLE resuhs'' are also tabulated in Table H 79 
of our 224 EROs are unmatched with any entry in the 
MDS database. 10 have more than 1 match within l". We 
include these objects and take the nearest match. YT03 
found that, for their sample, the MLE morphologies gave 
a mix of 50% disks and 37% bulges. Using the same scheme 
for our K ^ 18.7 subsample, we find relative fractions of 
32 ± 7% disks and 20 ± 6% bulges with 48% unclassified. 17 
of this subsample of 50 are not listed in the MDS catalogue 
jGriffiths et alJll994l) . Ignoring for the moment the objects 
without classifications, the MDS MLE classifications for the 

^ The Medium Deep Survey catalog is based on observations with 
the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA con- 
tract NAS5-26555. The Medium-Deep Survey analysis was funded 
by the HST WFPC2 Team and STScI grants G02684, G06951, 
G07536, and G08384. 



Table 3. Summary of the coarse morphological classifications for 
two of our ERO samples: the full sample and a subsample cut 
at K = 18.7, the median limiting magnitude of YT03. Two clas- 
sification schemes are considered: our visual morphologies and 
automated maximum likelihood estimator (MLE) fits. ^ - Galax- 
ies in our visual classification scheme which do not fit into either 
the disk or spheroid bins, and galaxies which are unclassified by 
the MDS MLE. See text for details. 



Sample 


Disks 


Spheroids 


Othert 


Total 


Visual 










K < 18.7 


27/54% 


12/24% 


11/22% 


50 


Full 


84/37% 


33/15% 


107/48% 


224 


MDS MLE 










K < 18.7 


16/32% 


10/20% 


24/48% 


50 


Full 


53/24% 


33/15% 


138/62% 


224 



YT03 sample find around 40% more disks than bulges; we 
find around 60% more disks than bulges. The general trends 
in these two surveys at bright {K ^ 18.7) are consistent: 
visual classifications identify approximately twice as many 
disks as bulges and the MDS MLE method identifies ap- 
proximately 60% more disks than bulges for those objects 
sufficiently bright to classify. 

We have shown via direct comparison that our visual 
morphologies are in reasonably good agreement with those 
of YT03, with a disagreement at the level of <, 20%. The rela- 
tive morphological mixes within both samples (YT03 versus 
our bright subsample) also appear consistent through both 
visual morphologies and automated bulge-|-disk decomposi- 
tion. Thus we can now use our full sample to look at evo- 
lution across the break in the ERO number counts into the 
K = 19-20 regime. 

3.3.6 Comparison with other morphological studies 

Our morpholog i cal fr actions contradict the findings of 
iMoriondo et alJ ll2000l) who found 15% irregular / disk-like, 
~50-80% elliptical-like. Nevertheless, as they themselves 
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note, their sample is drawn from a heterogeneous collec- 
tion of a rchival HST pointi ngs and may not be repre- 
sentative. iRoche et alJ JiflQ^) Eilso find a high fraction of 
bulge-dominated EROs using ground-based K-ha.nA data 
(39 ± 11% bulge, 25 ± 9% exponential disk profile systems) 
from a subsample of 32 of their EROs with high resolu- 
tion UFTI data taken in good seeing. Unfortunately, ground 
based NIR data is not ideal for morphological studies, as 
even in the best seeing conditions, the very nature of the 
NIR observing strategy means that multiple, dithered im- 
ages must be registered and combined, with the added com- 
plexity that the seeing may be changing between exposures. 
Correct propagation of all such uncertainties through the 
surface brigh tness profile fitting procedure is non-trivial. 
Furthermore. iRoche et ahl ll2002ll report a possible z~l clus- 
ter of EROs comprising five bulge systems and two disks, 
thus infiating the early-type galaxy fraction. 

Our morphological fractions are also in broad agree- 
ment with thos e from an H ST imaging survey of lensed 
EROs given by ISmith et all ll2002ah (18% compact, 50% 
irregular/disk-like and 32% unclassified) using an (R — K) ^ 
5.3 selection criterion. Although, we caution that (R — K) ^ 
5.3 selection may find difi^erent systems from (/ — K) ^ 4.0 
selection (see t|4.71 . 



Bright (K<19,12) 

Faint (K>19.38) 



Figure 5. A comparison of the cumulative colour distributions 
of the 100 brightest (solid line) and 100 faintest (broken line) 
EROs in our sample. Error bars are based on Poisson statistics 
in each bin. The fainter sample is skewed toward redder colours, 
indicating that the colour distribution evolves, as galaxies pass 
over the A" ~ 19 break in the number counts. 



4.2 Evolution of ERG morphologies 

We now investigate whether the changes in the colours of 
EROs across the break in their number counts are mirrored 
by a change in their morphological mix. 



4 RESULTS & DISCUSSION 

4.1 Evolution of ERG colour distribution 

To search for changes in the ERO population across the 
break in the number counts, we will take subsamples of the 
100 brightest and 100 faintest EROs in the survey, to main- 
tain equivalent uncertainties in both bins. This corresponds 
to A' < 19.16 and K > 19.38 for the bright and faint samples 
respectively. This is a good choice as, although the position 
of the break is not precisely defined, it seems to occur around 
K ^ 19.0-19.5. 

The cumulative colour distributions of the ERO sample 
divided into bright and faint A'-band bins are examined in 
Fig. |S] The distributions clearly appear different such that 
the fainter sample is skewed toward redder colours (with a 
median colour of (/ - AT) = 4.39 ± 0.07 and (/ - A) = 
4.95 ± 0.09 for the bright and faint samples, respectively). 
A two-sided KS test reveals that the probability of the faint 
colour distribution being drawn from the same population as 
the bright sample is ~ 3 x 10"*. In this comparison we have 
also chosen to conservatively treat lower limits in (/8i4 — K) 
as detections. Since all but two of the non-detections occur 
in the faint bin (see Fig. 0, objects in the faint sample can 
only become redder than the values tested. This strengthens 
the conclusion of a reddening of the fainter EROs. If we 
completely excise the break region and repeat the test using 
K <19.0 (81 objects) and K >19.5 (86 objects), the trend 
is unchanged, and the significance only marginally reduced 
to ~6xl0-*. 

Thus, the colours of EROs evolve as they pass over the 
break in the cumulative counts, such that a greater fraction 
of the faint EROs are redder than their bright counterparts. 



4-2.1 Morphological number counts of EROs 

We give the distribution of the different morphological sub- 
classes of EROs on the colour- magintude plane in Fig. |Sf . 
The main point to note from this plot is that the bright- 
est EROs, which are responsible for the sharp rise in the 
number counts at AT <: 19, are dominated by galaxies with 
disk-like morphologies. Indeed at AT <: 18, Z 80% of the 
sample are disk systems. To better illustrate the varia- 
tion in the morphological mix with magnitude we examine 
the morphologically-classified number counts in Fig. |7| At 
A" ~ 18 bulge systems start to appear, with an even steeper 
slope than the bright disk population. Merging systems also 
appear at A" ~ 18 and seem to follow the bulge population. 
To quantify these comparisons, we fit power laws to the mor- 
phological subsamples for K ^ 19.0, where the visual clas- 
sifications are reasonably complete (~ 90%), which yields 
gradients of adisk = 0.77 ± 0.11; abuige = 1.17 ± 0.27. For 
the unclassified (class 0) population: Qundass = 0.57 ± 0.19 
for K s; 19.0, and Oundass = 0.76 ± 0.40 for 19 < AT < 20. 
The remaining classes contain too few objects and/or are 
not well fit by a power law in this magnitude range. 

^ One particularly bright, red source which is unclassified is 
ER0198 with (I - K) = 6.05 ± 0.09 and K = 17.87 ± 0.03. 
This source lies on the edge of a halo from a bright star. It ap- 
pears to be a genuine detection, and has an F814W counterpart, 
but its proximity to the star halo may bias the magnitudes mea- 
sured, as the local sky background possesses a large gradient. 
However, multicolour photometry in additional bands yields a 
plausible photometric redshift of z ~1.7 for this object (see i|4.5l 
and it seems to just be a high redshift, extreme ERO. Note that 
similar objects are seen in the YT03 data. 
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Figure 7. Left panel shows the morphological number counts for our data and the double power law fit to the total ERO counts. We do 
not attempt to correct these counts for incompleteness as a function of morphological type, as this correction is negligible compared with 
the Poisson uncertainties we assume throughout. Note that nearly all the brightest EROs are disks, and bulges and mergers both appear 
at K ^ 18. Right panel shows morphological number counts for all K-sclectcd galaxies using MDS MLE classifications from F814W data. 
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Figure 6. The colour-magnitude diagram for the EROs labelled 
by their visual morphologies. Note the dominance of disk-like 
morphologies in the brightest EROs in our survey, with bulge 
systems only starting to appear a.t K > 18. 

Summarising these results in the scheme of Table |21 to 
K = 20, our visual scheme gives 37 ± 3% disk-like EROs 
and 15 ± 2% bulge-like, with 48% not fitting easily into 
either of these simple classes (19% of which are unclassi- 
fied). Again, we see around twice as many disks as bulges, 
but now the unclassified-l-other fraction is almost twice as 
high as for the bright subsample. Similarly if we use the 
MLE classifications, we find that the ratio of bulge to disk 
galaxies is consistent within the uncertainties with that mea- 
sured for the bright subsample, but the relative fraction of 
unclassifiable-|-other galaxies has doubled. 

4-2.2 Comparison with morphological number counts of 
all K-selected galaxies 

We also fit power laws to the /f-selected morphological 
counts for all galaxies using MDS MLE morphologies from 
F814W data (Fig.[7|l. We find slopes of q = 0.30 ± 0.02 for 



the disks and a = 0.35 ± 0.02 for the bulges, in the range 
18.0 ^ K ^ 20.0. The unclassified population follows the 
bulge counts closely. The counts for these sub-populations 
are substantially shallower than their counterparts in the 
ERO sample. Our preferred interpretation of the steepness 
of the ERO counts is that the EROs are a high redshift 
(z~ 1, inferred from a passively evolved M* for the bright 
EROs of K~18) subsample of the total is'-selected galaxy 
population. We return to the question of redshifts in t|4.5l 

4.2.3 The nature of the bright disk EROs 

To better understand the nature of the bright, ERO disk 
galaxies we first measure the average semi-major axes of 
all the well-detected disk systems (as measured by SEx- 
tractor's Ajmage parameter). We find a mean of <A>= 
(0.16 ± 0.01)"; while the brighter disk EROs (K< 18) are 
also systematically bigger: < A >= (0.32 ± 0.03)". Thus 
these bright (and relatively blue) disk EROs also appear to 
be larger than the average morphologically-classified disks 
in the ERO population. Their elongations range from 1.12 
to 1.61 with a median of 1.36, suggesting that these systems 
are only moderately inclined 45° to line of sight). If these 
EROs are at z~l-2, then the median size of the bright disk 
sample implies a physical half light radius* of ^ 3kpc, com- 
pared to ~ 2kpc for the full sample of disk EROs. These sizes 
are comparable to those seen for similarly luminous dis k 
galaxies in the local Universe (e.g. iBurgarella et all 1200 iT) . 
Thus, although these disk EROs are among the brightest 
and largest galaxies in the sample, they are neither too 
bright nor too large to be conclusively ruled out as ly- 
ing at z>l. Indeed, recent work has uncovered large disks 
with half light ra dii of 5.0-7.5kpc) at redshifts z=l. 4-3.0 
Labbe et ailbOOSTl . 

YT03 claimed that around 40% of their disk-like EROs 
(or one third of all their EROs) showed unusually large, edge 
on disks and suggest that these systems are actually lower 

* for all half light radii we use the semi-major axis 
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redshift contaminants to the 1 ERO population. We find 
only one very large edge on disk ERO of comparable size 
to the one illustrated in YTOS's fig. 8 (EROOll with a half 
light radius of 0.45", Fig. IIH . Thus we suggest that this 
contamination is slight for the 20 ERO population. 

4-2.4 ERO sub-populations 

Next, we examine the colours and apparent magnitudes of 
the various morphological sub-populations of EROs in Ta- 
ble |1] We consider the disk-like and bulge-like classes, the 
unclassified population, the full sample and the amorphous 
object s (these should be clos est to the irregular EROs stud- 
ied bv lMoriondo et al]|2000h . As would be expected, we see 
that the median if-magnitude of the unclassified population 
is fainter (by >0.5 mag) than any other class. The median 
colour of the faint sample for the 'all' class is seen to be 
redder by ~ 0.3 mags than the faint sample. This is the 
colour evolution reported in H4.ll The median colours of 
both the bright and f aint disks and bulges are consistent. 
iMoriondo et alj i2000l) found that irregular EROs were red- 
der than their bulge EROs. We see that our amorphous sys- 
tems are indeed redder than our bulges, particularly in the 
bright subsample. We also note that in contrast to the full 
sample, in this case it is the bright amorphous subsample 
which is redder than its faint counterpart. 

Finally, we briefly examine the environment of the 
EROs, as defined by the number oi K ^20 galaxies within 1 
arcmin of each ERO. This radius corresponds to a physical 
size of ~0.5Mpc at z~l. We failed to find any clear trend, 
with all morphological types inhabiting similar density en- 
vironments. 

4.3 Quantitative evolution of ERO morphologies 

In order to quantify the form of this evolution, we artificially 
faded our bright subsample of EROs to the magnitudes of 
our faint subsample and statistically compared properties of 
the artificially faded population with the true faint popula- 
tion. The first quantitative measure we examine is concen- 
tration, C (ET^I. 

4-3.1 Artificially fading the bright ERO population 

To artificially dim the EROs, each of the 100 A'-brightest 
EROs has an i814 magnitude selected randomly from the 
observed distribution of the 100 A'-faintest EROs. Each 
bright ERO was then extracted from the HST image, 
dimmed to its corresponding faint magnitude (keeping a 
fixed angular size) and reinserted into a random blank re- 
gion within the original F814W image. Its concentration was 
then remeasured. By repeating this for each of the 100 bright 
EROs, an artificial faint sample is created which can be sta- 
tistically compared with the true faint sample. By directly 
fading a bright subsample of our objects, all the uncertain- 
ties associated with measuring properties of the faint sam- 
ple are accounted for. These simulations were repeated 100 
times and the maximum probability of the two samples be- 
ing drawn from the same parent population was found to be 
0.6%, with the mean probability of the 100 simulations being 
0.1%. Thus, we detect evolution in the concentration indices 



for EROs across the break in the ERO number counts. Note 
that the angular size change over the expected redshift range 
of the objects is unlikely to be responsible: the angular size 
of an object only changes by ~ 5% between a redshift of 1 
and 2. Again, this result still holds if we change the bright 
and faint samples limits to K<19.0 and K>19.5, although 
the maximum probability is increased to 14% and the mean 
is 3%. 

4-3-2 Evolution in ^ — C space 

The surface brightness - concentration index (/i — C) classi- 
fication plane has been used to quantify galaxy morphology, 
with bulge-dominated galaxies typically residing at higher 
concentration and greater s urface brightnesses th an disk- 
dominated objects (see fig. l. lAbraham et allll994^ . 

Fig-IHshows the distributions for our 100 brightest and 
100 faintest EROs. Representative error bars in C deter- 
mined from simulations are indicated at the top of each 
panel and for clarity the errors in surface brightness are 
shown only for the bulges. The spread of sources in the 
bright ERO sample is greater than the estimated errors in 
both quantities, indicating that these quantitative observ- 
ables are measuring a real variation in the morphological 
properties within the EROs population, although the dis- 
tributions for the visually classified bulges and disks over- 
lap significantly. We can identify a region in this plane (de- 
scribed below) towards higher concentration and higher sur- 
face brightnesses which contains the bulk of the visually 
identified bulge-like population. We take the locations of our 
visually classified bulges and disks in — C space as further 
support for the reliability of our visual classifications. In ad- 
dition, we note that peculiar (i.e. merger and amorphous) 
systems do not occupy specific areas of the diagram and are 
scattered throughout all the occupied range. 

To test the reliability of the apparent evolution on the 
H ~ C plane, we concentrate on the bulge population and 
perform a simple comparison of the expected numbers of 
faint, bulge-classified EROs. In the left panel of Fig. |S]we 
identify a rectangular region which contains a large fraction 
of the visually identified bulges. We define the completeness 
and purity of this box as the fraction of bulges out of the 
total number of bulges (in this sample of 100) located in 
this box, and the fraction of bulges out of the number of 
objects within this box, respectively. For the 35 objects in 
the box in the bright subsample we find a completeness of 
0.76±0.11 and a purity of 0.46±0.19. In the right panel, we 
move the box by the mean fading vector determined from 
the simulations, and recalculate these values for the faintest 
100 EROs, finding a completeness of 0.42 ± 0.16 and purity 
of 0.36 ± 0.19 for the 53 objects now in the box. The box 
selected is quite large and in the bright sample is contami- 
nated by at least as many non-bulges as bulges. However, the 
translation of this box and the location of visually classified 
bulges in the faint sample shows that the bulge population 
behaves as expected, moving along the simulated fading vec- 
tor, and a comparable fraction of the faint bulge population 
is still located within this region. We thus conclude that the 
changing morphological mix between the bright and faint 
ERO population is not a result of visual misclassification. 
The contours in the right panel denote the density of ar- 
tificially faded EROs, the bulk of which lie toward lower 
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Table 4. Average properties of ERO sub-populations. We give the median colours, and median V^Tot magnitude for various morphological 
sub-populations. Amorphous objects appear significantly redder than the bulge population (particularly in the bright subsample) as noted 
by Moriondo et al. (2000). Note also that the faint subsample for all morphological types is redder than the bright subsample. This trend 
is reversed when just considering the amorphous systems. Errors are uncertainties on the median from bootstrap resampling. 



Sample 


Disks 


Bulges 


Amorphous 


Compact 


Merger 


Unclassified 


All 


Median(/ - K) 

K < 19.0 

19.0 ^ i^: < 20.0 


4.30 ±0.05 
4.41 ±0.10 


4.40 ±0.17 
4.40 ±0.12 


5.18 ±0.31 
4.73 ±0.20 


4.19 ±0.08 
4.29 ±0.17 


4.58 ±0.22 
4.57 ±0.13 


6.68 ±0.41 
5.54 ±0.10 


4.38 ±0.05 
4.71 ±0.08 


#/% 

K ^ 19.0 

19.0 <^ 20.0 


43/53.75% 
41/28.47% 


15/18.75% 
18/12.50% 


2/2.50% 
16/11.11% 


4/5.00% 
28/19.44% 


9/11.25% 
5/3.47% 


7/8.75% 
36/25.00% 


80/100% 
144/100% 


Media.n{KTot) 


18.97 ±0.12 


19.04 ±0.18 


19.56 ±0.08 


19.54 ±0.12 


18.91 ±0.18 


19.63 ±0.10 


19.27 ±0.06 




Figure 8. The distribution of mean i814 surface brightness (/i) 
versus concentration index (C) for the 100 brightest (left panel) 
and 100 faintest (right panel) EROs in our survey. Symbols de- 
note visual morphologies and are the same as for Fig. El The 
uncertainty in /i is just shown on bulge classes for clarity and 
the typical errors in C estimated from simulations are shown at 
top of each plot. The box in the left panel is the region selected 
to contain a large fraction of visually identified bulges (^ 70%). 
The arrow in the right panel indicates the translation in ii—C ob- 
served when the bright ERO population is faded to match the 
faint population (as described in text) and the dashed box shows 
the movement of the translated classification box. Contours in 
the right panel show the density of artificially faded galaxies (of 
all morphological types). It appears that the distribution of the 
simulated faded population extends to lower C than is seen in the 
real faint sample, suggesting that the low-C galaxies seen in the 
bright sample are absent in the fainter population. 



values of C than seen in the real faint sample (as measured 
in S lOTl . 

Hence, it would appear that morphological changes are 
associated with a decline in the number of low concentration 
objects. Galaxies in this region of the fi — C plane appear 
to be disks in the bright sample (a nd would be expecte d 
to be disks in the local universe, e.g. lAbraham et alJIlQQj) . 
Thus we interpret the break in the ERO counts as a real 
decrease in the number of disk-like galaxies in the population 
at K ^ 19. 



Figure 9. The (/814 ~ K) versus (J — K) colour-colour diagram 
for subsample of our EROs with J-band photometry. EROs are 
labelled with their visual classifications. The solid line shows the 
proposed division between photometrically classified passive and 
star-forming galaxies (Pozzetti & Mannucci 2000). This division is 
broadly supported by our morphological classifications, however a 
number of EROs with disk-like morphologies appear in the lower 
left corner and hence fall in the elliptical category. As we show 
in H4.4I a more detailed modelling of the SEDs of these galaxies 
suggests that they are star-forming galaxies which this simple 
photometric classification scheme is failing to correctly identify. 

4.4 Photometric classification of EROs 

Another possibility for distinguishing evolution in the 
various subclasses o f ERO s is multicolour photometry. 
iPozzetti fc Mannucci i200Cf) proposed a photometric clas- 
sification scheme for separating dusty star-forming SEDs 
from passive systems in the I/J/K colour plane. Us- 
ing the subsample of our survey with J-band photome- 
try we illustrate this in Fig . |5| The solid line denotes the 
iPozzetti fc Mannucci feOO d) division between the bluer, pas- 
sive and redder, star-forming systems. We can compare this 
classification scheme with the distribution of our visual mor- 
phological classifications to test the extent of the overlap 
between these two approaches. Here we start from the naive 
assumption that all morphologically-classified bulge systems 
will have passive SEDs, while any galaxy with a disk com- 
ponent is probably an ERO by virtue of a strong dust com- 
ponent produced by ongoing star formation activity. 
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51 EROs lie in the passive region of Fig. |5| and 21 in 
the star-forming region. Firstly we quantify the fractions 
of photometrically classified objects in each morphological 
class. Of the 30 disks, 21 (70%) are classified as having star- 
forming SEDs, with only 9 (30%) falling into the passive 
region. Equally, 10 (91%) of the 11 bulges are classed as pas- 
sive and 11 (78.6%) of the 14 compact sources have colours 
consistent with passive SEDs. For the remaining classes, 2 
(66%) of the 3 mergers; 3 (50%) of the 6 amorphous and 6 
(60%) of the 10 unclassified objects have SEDs classed as 
star- forming. 

Thus the photometric classifications crudely support 
our prejudices based on the morphologies of the EROs: 
most disk systems have star-forming SEDs and most bulge 
systems have passive colours; indeed we find only one 
morphologically-classified bulge lies in the star-forming re- 
gion. The main discrepancy comes from a number of disks 
falling into the passive region of the plot, although these are 
mostly relatively blue in (/8i4 — K), with (/8i4 — K) ^ 4.5. 
Closer inspection of the morphologies of these EROs does 
not show strong evidence of them being early-type disk 
galaxies (which could be expected to show passive SEDs). 

Another clear trend from this comparison is that a simi- 
lar fraction of compact sources are classed as passive to that 
seen for the bulge population. This is consistent with the 
idea that the compact class comprises bulge systems which 
are too faint to be unambiguously classified as such. The 
amorphous, merging and unclassified samples show much 
higher proportions of sources with star-forming SEDs, sug- 
gesting that these indeed represent morphologically more 
complex systems, rather than misclassified faint bulges. 

Hence using our visual morpholog ies we find good agree- 
ment with the photometric scheme of lPozzetti fc Mannuccil 
l|2000l) for the higher signal to noise bulge-classified galaxies 
and less agreement for the disk-classified ERO sample. We 
use this to infer the probable morphological mix of the lower 
signal to noise amorphous and compact EROs, associating 
the amorphous systems with star-forming and the compact 
systems as mostly passive. 

S02 used the equivalent R/ J /K colour-colour diagram 
to analyse the colours of radio-detected and undetected 
EROs. They found a similar level of agreement, in that 
the EROs expected on the basis of their radio emission to 
be star-forming mainly fell into the star-formi ng region of 
the dia gram, with some cross-contamination. ISmith et alJ 
(|20023) also used this colour-colour diagram to examine the 
location of their HST imaged EROs. They found that all 
their objects with compact morphologies lay within the pas- 
sive region. The size of their J-band photometric errors pre- 
vented them from drawing conclusions from their irregular 
objects. An alternative but similar colour selection crite- 
rion in (I — J)-(J — K) colour-colour space was used by 
IWehner et alJ ll2002h to study the colours of EROs with 
sub-mm observations. They found that objects in the star- 
forming region of the diagram were indeed responsible for 
the bulk of the 850 ^.m emission detected from their coad- 
ded observat i ons, in dicative of dusty star-formation activity. 
iRoche et alj (|20o3) adopted a similar scheme using galaxy 
tracks plotted in 1/ J/H/K colour space to attempt to sepa- 
rate the classes. They noted, as with these other works, that 
a large fraction of the EROs lie close to the dividing line be- 
tween classes. This could be symptomatic of the size of the 



photometric errors, or indicate that most of this class of 
object do indeed possess SEDs intermediate between dusty 
star-forming and passive stellar populations. 

4.5 Photometric Redshifts 

22 of our ERO sample have good detections in four pho- 
tometric bands; Reoe and /8i4 from the MDS observations 
and J and K from our INGRID survey. The MDS F606W 
observations of these fields were retrieved from the ST-ECF 
archive and reduced in a similar manner to the F814W im- 
ages ( t|2.2ll . For these galaxies it is possible to investigate 
the properties of their SEDs in more detail. However, the 
apparent colours are affected both by their intrinsic SED 
and their redshift and we are therefore required to fit for 
both of these variables. H ence, we have used a photomet- 
ric redshift code, hyperz iBolzonella et al.ll2000h . to study 
these objects. We adopted a procedure similar to S02. We 
attempt to fit, in turn, a dusty star-forming and an evolved 
SED to each object. We use a single star-formation history 
with an e-folding time-scale of r =lGyr, a Miller-Scalo IMF 
and solar metallicity. The dusty SED is allowed reddening 
values in t he range Ay = 1 — 6 (consistent with values es- 
timated bv lCimatti et al]|2002h . and the evolved SED must 
have a reddening value Ay ^ 0.2. Redshifts in the range 
to 4 are considered. We list the photometric data and the 
derived properties for this sample in TableO In this analysis 
we retain the F606W passband, rather than introducing any 
uncertainty by transforming it to another, more common, 
filter system. We consider only the EROs with detections 
in all four passbands, rather than including limits as well, 
as using three or fewer detections leads to very poorly con- 
strained SEDs. If this selection introduces any bias, it will 
be that requiring detections in the optical passbands pref- 
erentially selects bluer galaxies. However, we are primarily 
interested in the relation between SED classifications and 
colour/morphological classifications, rather than construct- 
ing a fair sample. 

From the 22 EROs, 14 could be fit using hyperz; the re- 
maining 8 are typically the reddest objects in (F606W— Tf). 
All but one have {F606W-K) ^ 7. The median redshift 
of our sample with photometric redshifts, which has a me- 
dian magnitude oi K = 19.43 ± 0.19, is z= 1.20 ± 0.22 
with al l but two E ROs predicted to lie in the range 0.7 
to 1.8. ICimatti et a l. (2002) found a median redshift of 
<z>— 1.1 ± 0.2 with galaxies in the range z = 0.7-1.4, 
for [R- K) > 5.0, K ^ 19.2 EROs, using optical spec- 
troscopy. Thus our samp l e spa ns the same redshift range 
as that of ICimatti et al.l (|2002|) . but with a higher upper 
limit, as might be expected from the fainter magnitude range 
probed here. 

For the fitted EROs we find 11/14 have disk-like mor- 
phologies, with a further 2/14 either amorphous or LSB 
and one unclassified (the high fraction of morphologically- 
classified EROs in this relatively faint A'-band subsample 
simply reflects our selection of EROs with strong optical 
detections) . 

The photometric classifications give a dusty star- 
forming SED as the best fit for all but two cases which favour 
evolved SEDs. Interestingly, half of these EROs (7/14) with 
dusty star-forming SEDs lie in the passive region of the 
(/814 - K)-{J - K) plane (Fig.EJ, with many of the EROs 
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occupying the (/8i4 — K) ^ 4.0 region in the lower left cor- 
ner of the diagram. This may reflect our selection, where by 
requiring the EROs be red in (7 — K) but sufficiently blue 
to be well detected in I^u and F606W, we select unusual 
systems, possessing strong UV upturns and possibly a mix- 
ture of young and evolved stellar populations. We note that 
of the blue disks in the passive region, 29% (2/7, or 2/9 if 
we include two borderline passive cases) are best fitted with 
evolved SEDs, suggesting that these are indeed passive disk 
systems. 

One of these EROs (#226, an amorphous object) is 
best-fitted as a low-luminosity, dusty system at low redshift 
(z^O.4), but this is only marginally favoured over an evolved 
SED at z=1.2. Overall, it seems unlikely that a large frac- 
tion of our EROs are at such low redshifts ( z'^0.4). More- 
over, s uch systems are not seen in the brighter lCimatti et alJ 
||2002) spectroscopic sample. 

4.6 High-redshift EROs 

One of the EROs in the photometric analysis in the pre- 
vious section has a much higher estimated redshift than 
the remainder of the sample. This ERO is #158 at Zphot ~ 
3.39lQ2g. This source is one of only two EROs in this sub- 
sam ple with (J — K) ^ 2.3. Van Dokkum et al. 2003 (see 
also lFranx et al.(l2003l) have recently shown through optical 
spectroscopy that selecting galaxies with (J — K) ^ 2.3 effi- 
ciently selects objects with a prominent optical break (either 
the 3625A Balmer break or the 4000A Call H-f K break) at 
z ^ 2. We can use all of our MDS fields with I / J/K photom- 
etry (81 arcmin^) to examine the nature of this class of ERO. 
Adopting this selection, we find 24 EROs with (J—K) ^ 2.3, 
or a surface density of (0.30 ±0.06) arcmin~^ (Poisson error) 
to K ^20. Van Dokkum et al. 2003 only give the surface den- 
sity at if = 21 {l.OQtoTe) 

and so we are unable to compare 
with their numbers, except to note that our surface density 
is lower to this brighter magnitude limit, and that there are 
very large field to field variations in this class of ERO (there 
are n o (J-K) ^ 2.3 EROs in the HDFS. Ivan Dokkum et "ail 
l2003l) . The morphological mix of these galaxies is 20 ± 9% 
disk-hke, 4 ± 4% bulge-like, 17 ± 7% amorphous, 8 ± 6% 
merger, 13 ± 7% compact and 38 ± 13% unclassifiable. The 
fraction of bulges relative to disks in the sample is lower 
than for the full {I — K) ^ 4.0 ERO sample, although just 
compatible within the 1-a Poisson errors. The fraction of 
peculiar (not disk-like or bulge-like) EROs is higher than 
for the full sample. The magnitudes of these (J — K) ^ 2.3 
selected EROs span the full range of the (/ — K) ERO sam- 
ple. It is interesting to note the lack of bulge-like systems 
in this perhaps higher r edshift sample, and also t he absence 
of passive spectra from Ivan Dokkum et al. ' (2003)'s sample. 
This could be an indication that at these redshifts (z^ 2) the 
formation epoch of early type galaxies is being approached 
(as suggested by complimentary studies of cluster ellipticals, 
e.g. [Stanford ct al 199&) . 

4.7 ERO sample differences in different 
photometric passbands 

A range of different photometric criteria have been used to 
identify 'Extremely Red Objects'. How do the different cri- 
teria affect the mix of objects selected? The ERO number 



counts for samples selected by a variety of definitions (e.g. 
(R — K) J5 5.3, (/ — H) ^ 3.0) are in broad agreement (e.g. 
fig. 7 of YT03), and in addition the K « 19-20 break in the 
counts is seen in sufficiently deep dat a selected by (Ign — K) , 
(I - H), and (,^702 - K) (this work:' McCarthv et aljliooil : 
ISmith et alJl2002al . respectively) . However, a direct compar- 
ison of the details of the samples such as the relative mor- 
phological mixes of EROs using different colour selection has 
yet to be addressed. 

The difference between the {R — K) and (/ — K) 
colour cuts can be examined for a subsample of our ERO 
sample with F606W data. F606W can be approximately 
transformed to the Cousins R passband using Reoe ~ 
F606 Wveoa - 0.37(F606W-F814W)ve£,a JMetcalfe et alJ 
120011) . This results in 27 EROs with measurable J?606 pho- 
tometry (Fig. not . Only three of the (78i4 — K) selected 
EROs are bluer in {Rsoe — K) than the typical selection 
criterion of (i?606 — K) ^ 5.3. One of these sources is very 
close to the selection boundary a nd would be found using 
the {R(,QQ—K) ^ 5.0 criterion ("e.g. lCimatti et all2002ll . and 
may even be found with the former criterion, given the un- 
certainty in the F606W to Tieoe transformation. The remain- 
ing two EROs have (7?606 — K) ~ 4.0 and would not. Thus, 
~90% of our EROs would also be selected in an (i?606 — K) 
ERO survey (in agreement with S02 who found a value of 
93%), suggesting that (78i4 — K) ^ 4.0 is a less stringent 
ERO requirement than (i?606 ~ K) ^ 5.3. Furthermore, for 
this admittedly small sample, all the galaxies dropping out 
of the ERO category when selected in [R — K) are disks, 
which supports the suggestion by YT03 that (I — K) selected 
ERO surveys may preferentially include disks. To quantita- 
tively compare the different morphological mixes found by 
the two techniques requires R and / band data and HST 
morphologies on the same regions of sky (e.g. Gilbank et al. 
in prep). 

We also performed visual classifications of the EROs in 
F606W to compare with the morphologies derived from the 
F814W images. We find for the 27 EROs with sufficient sig- 
nal for F606W photometry: 8 of these are too faint to mor- 
phologically classify; 10 of these are disk-like in F606W and 
were also disk-like in F814W; 2 galaxies were compact and 
disk- like in F814W and were simply compact in F606W; 5 
galaxies appear LSB and disturbed in F606W whereas these 
systems were just classed as disk-like in F814W; one galaxy 
appears amorphous in F606W and disk-like in F814W. One 
ERO classed as compact and disk-like in F606W was classed 
as compact and symmetrical in F814W (nominally a bulge 
class) . This represents the only discrepancy between our no- 
tional bulge and disk classes, but the difference is slight as 
in both passbands the object is compact and symmetric. 
Thus, the agreement between the F814W and F606W vi- 
sual morphologies is very good. There is a tendency for the 
disks in the latter to exhibit more disturbed morphologies, 
as would be expected as the F606W probes further into the 
extreme ultraviolet, in each EROs' rest-frame. This agree- 
ment between passbands is reassuring, as F606W at z~ 1.0 
corresponds to the same rest-frame wavelength as F814W 
at z~ 1.5 - a plausible redshift spread for the ERO sam- 
ple. For this small subsample, we only find one (marginal) 
bulge - disk misclassification, whereas we assume Poisson 
uncertainties on all our measured fractions throughout this 
work. 
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Table 5. Photometric properties of EROs with F606W, /814, J and K data. In the SED fitting procedure, the raw F606W magnitudes 
have been used directly instead of attempting to convert to -Reoe ■ " — best fitting SED type: D — dusty; E - evolved. ^ - PMOO shows 
the classification based on the Pozzetti &; Mannucci scheme. P — passive, S - star-forming. Square brackets indicate systems close to the 
dividing line, whose errors allow them into the other class. - the 99% confidence interval for the photometric redshift. - morphology 
refers to visual classification. 



ID (F606W-ft:) (/814 - K) {J - K) K SED" PMOO'' Zp^^t Zp^ot Morph'' Comments 

(99%)'^ (visual morphology) 



001 


6 


93± 





14 


4 


60± 





13 


2 


05± 





22 


19 


69 


D 


002 


6 


62± 





13 


4 


33± 





13 


1 


37± 





16 


19 


63 


D 


158 


5 


78± 





06 


4 


13± 





05 


2 


84± 





06 


18 


11 


D 


167 


6 


04± 





14 


4 


13± 





13 


2 


29± 





17 


19 


73 


D 


197 


5 


66± 





09 


4 


20± 





09 


1 


42± 





14 


19 


87 


D 


198 


8 


78± 





23 


6 


05± 





09 


2 


53± 





15 


17 


87 


D 


204 


7 


27± 





17 


5 


25± 





11 


2 


17± 





23 


19 


58 


D 


206 


7 


27± 





05 


4 


30± 





03 


2 


08± 





09 


18 


78 


D 


209 


8 


42± 





10 


5 


25± 





04 


1 


95± 





07 


18 


48 


D 


224 


5 


77± 





07 


4 


18± 





06 


1 


17± 





10 


19 


43 


D 


225 


6 


18± 





05 


4 


17± 





04 


1 


27± 





07 


18 


53 


E 


226 


5 


84± 





08 


4 


19± 





07 


1 


40± 





15 


19 


65 


D 


227 


5 


97± 





07 


4 


72± 





05 


1 


46± 





11 


18 


63 


D 


229 


6 


09± 





06 


4 


07± 





05 


1 


29± 





09 


19 


14 


E 



[S] 


1 


20 


0.28- 


-1.49 


2 




P 


1 


18 


0.93- 


-1.51 


2 


early-type? 


s 


3 


86 


3.20- 


-4.00 


6 


edge-on late-type disk 


p 





71 


0.15- 


-1.26 


6 


small edge-on disk 1"/ 
late-type 


s 


1 


49 


1.26- 


-1.79 


2 




[S] 


1 


71 


0.83- 


-1.93 







[P] 


1 


00 


0.95 


-1.9 


2 




s 





89 


0.58- 


-1.18 


6 


edge-on Scd/Sdm 2.5" 


[P] 


1 


49 


1.44- 


-1.66 


2 


(near chip edge) 


p 


1 


49 


1.34- 


-1.61 


6 


edge-on disk with tidal 
tail/cmp 2" 


p 


1 


09 


1.00- 


-1.10 


6 


face-on asymm disk 
1.5" 


p 





10 


0.00- 


-0.42 


9 




p 


1 


79 


1.51- 


-2.00 


5 


face-on disk or clumpy 
LSB 


p 


1 


06 


0.97- 


-1.09 


6 


late-type slightly in- 



clined disk 2' 



8 
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: + ^ 


^ : 

^ •Bulge : 
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: 
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Figure 10. The (Reoe — K)-{If,\i — K) colour-colour diagram 
for our ERO sample (selected by (/814 — K) ^ 4.0) with available 
F606W photometry. The dashed line shows the typical (iJsoe ~ 
K) ^ 5.3 selection bound for R-selected EROs, and the redder 
(^^606 ~ K) ^ 6.0 limit discussed in the text. Points are labelled 
according to their visual morphologies in broad classes - disk- 
like, bulge-like, or unclassified. Vertical error bars represent only 
the error in the F606W photometry and do not account for the 
uncertainty in the transformation to R^oe- This suggests that 
the (/814 ~ K) ^ 4.0 selection criterion is less stringent than 
{Reoa ~K)^ 5.3. 

4.8 Composition of the ERO population 

Our visual and quantitative analysis of the morphological 
mix of our ERO catalogue suggests that the ERO popu- 
lation at if ^ 19 is dominated by galaxies with disk-like 
morphologies (54% of the total sample), with EROs with 
bulge-like or compact morphologies making up around a fur- 



ther quarter and a modest fraction ('-^ 15%) of EROs with 
disturbed, merging or amorphous morphologies. As we have 
shown, this dist ribution is similar to that found by previous 
studies (YT03: ISmith et al.ll2002af) . In contrast, our faint 
ERO subsample, with K = 19-20, has a lower fraction of 
disk-like EROs, down to less than a third, with the fractions 
of the other classes effectively unchanged and the decline in 
the disk-like EROs being made up for by a corresponding 
increase in the fraction of those with amorphous, merger 
or unclassified morphologies. Our quantitative tests of this 
evolution confirm that the apparent decline in the disk-like 
fraction is real and is not due to our inability to classify faint 
examples of this population. 

We can now ask what limits we can place on the pro- 
portion of the ERO population arising from galaxies which 
are red due to dust-obscuration and what fraction comes 
from evolved, passive systems. We use our visual classifica- 
tions to divide our sample into star-forming or passive EROs 
by taking extreme values in the following way. We assume 
that the minimum fraction of star-forming EROs is given by 
all the disk-like galaxies (although note that passive disks 
are also hkely, S 14.51 Smith ct al. 2002b), and the maximum 
star-forming fraction is given by all the galaxies which could 
plausibly harbour star-formation (including the disks, amor- 
phous, mergers and all the unclassified objects). In the same 
way we associate all bulge-like EROs with passive systems 
and add all compact and unclassified sources to this to gen- 
erate an upper limit. We thus find 54-77% star-forming and 
19-31% passive EROs in the bright subsample {K 19.0), 
and 29-71% star-forming and 13-58% passive in the faint 
sample (19.0 ^ K ^ 20.0). For the full sample these num- 
bers are 38-73% and 15-48% for star-forming and passive 
EROs, respectively. S02 (based on an ERO sample selected 
using {R — K) ^ 5.3) used deep radio data to estimate the 
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fraction of star forming EROs, finding a value of ~ 45%, 
to K = 20.5. This is compatible with our estimate based 
on morphologies. Furthermore, they used photometric SED 
fitting to assess the fraction of dusty galaxies, placed a firm 
lower limit of ^ 30%, and estimated that possibly 60 ± 15% 
of this ERO population could be dusty star-forming galax- 
ies. Again, this does not contradict our findings. 

In the discussion above, it is the increasing fraction of 
morphologically unclassified EROs in the faintest samples 
which produces the greatest uncertainty in the evolutionary 
trends we derive. These galaxies make up a quarter of the 
K — 19-20 sample and appear to also be responsible for 
the reddening of the whole ERO population at K ^ 19, 
especially when balanced against the declining fraction of 
bluer, disk-like EROs (Table |1J. It is tempting to associate 
this population with the counterparts to the dusty, active 
source s selected in the sub- mm waveband (e.g. ISmail et alJ 
l2002al: IWehner et al.ll2002^ . However, the confirmation of 
this will have to await complimentary multiwavelength and 
morphological studies, with the latter either benefiting from 
the higher red sensitivity of the new ACS camera on-board 
HST and the revival of the NICMOS near-infrared imaging 
capability of HST. 

4.9 Cosmological significance of EROs 

Assuming the median photometric redshift (z 1.5) for our 
sample, the break around K = 19 corresponds to a lumi- 
nosity of Mv ~ —19.7. For the star-forming population, we 
correct this value for reddening to g ive My ~ —22.1 (as- 
suming Av ~ 2.4. ICimatti et al.ll2002r . This unobscured lu- 
minosity corresponds approximately to an Ly galaxy today. 
Conservatively assuming no further star-formation activity 
in these galaxies between z ~ 1.5 and the present day (and 
no substantial merging), these galaxies would actually cor- 
respond to sub-L^ galaxies a,t z — 0. Similarly, if we adopt 
passive evolution models for the passive/bulge population, 
then K = 19 a,t z ^ 1.5 corresponds to an Ly elliptical 
today. 

Taking the range of our photometric redshifts (0.7 ^ 
z ^ 1.8, ignoring the two outliers) we derive a comoving 
volume for our survey of ~ 2 x lO^Mpc"'^. This gives a 
space density for our fuU ERO sample of ~ 1 x 10"^ Mpc'^ 
This is only a crude estimate since we do not know the true 
redshift distributions of the galaxies, and furthermore, the 
redshift distributions of the dusty star-forming and passive 
EROs are unlikely to be the same. Thus we use only an order 
of magnitude estimate. This number density agrees to this 
level with that derived locally from SDSS-I-2MASS data for 
> L* galaxies ( Bell ct al. 2002^. Thus, the ERO population 
could plausibly account for a sizable fraction (^ 10%, and 
potentially all) of the stars seen in luminous (~ L*) galaxies 
locally. 

A more detailed comparison of the ERO counts with a 
variety of evolution models was undertaken by _S.r nith e t al. 
(|20023), who compared their ERO counts with Pure Lu- 
minosity Evolution (PLE) and semi-analytic hierarchical 
galaxy formation models. The only model to correctly repro- 
duce the cumulative ERO number counts was one in which 
the ERO population was modelled as a single population of 
passive galaxies. This model now clearly overestimates the 
number of passive EROs as we have been shown that a siz- 



able fraction of the EROs are not si mple passive systems . 
The more realistic models examined bv lSmith et al.l ll2002a^ 
- PLE with a mixture of galaxy types and semi-analytic 
galaxy formation models - both under-predicted the num- 
ber of EROs by around an order of magnitude. A compari- 
son with these models is beyond the scope of this paper and 
will be presented in future work. However, we note briefly 
that more recent semi-analytic models, which inc lude bursts 
of sta r formation at high redshift absent in the ICole et all 
(l200Cfl reference model, now reproduce the number counts 
for [I - K) ^ 4.0 EROs (C. Baugh, priv. comm.). Unfor- 
tunately, the numbers of objects in redder subsamples (e.g. 
(/ — K) ^ 5.0) are still underestimated by around an or- 
der of magnitude. The requirement of these starbursts at 
higher redshift (which is also indicated by the presence of 
luminous, passive galaxies at z ~ 1) is consistent with the 
interpretation of SCUBA sources as protoelliptical galaxies 
in the process of formation (e.g. ISmail et al.ll2002all . 



5 CONCLUSIONS 

We have compiled a new catalogue of if-selected Extremely 
Red Objects in regions with deep HST /gi4-band obser- 
vations. We apply a number of tests to examine the na- 
ture of the ERO population and investigate the origin of 
the turnover in the ERO number counts around K ~ 19. 
The main aim of our work is to study the morphologies 
of our EROs using the exquisite resolution of the HST. 
However, for a subsample of our dataset with multicolour 
data we can also utilise SED-fitting as a classification tool. 
Hence, we were able to examine the commonly applied hy- 
pothesis that the passive population defined by multicolour 
photometry (or equally spectral line diagnostics) is associ- 
ated with bulge-dominated morphologies; and that EROs 
with the colours of dusty star-forming galaxies have disk- 
likc/irregular morphologies. 

1. We find a surface density of (/gi4 — K) 4.0 EROs 
(1.14±0.08) arcmin"^ aX K = 20.0. The number count slope 
fiattens from a = 0.88±0.09 for K ^ 19.0 to a = 0.42±0.19 
for 19.0 ^ if ^ 20.0, in good agreement with other surveys 
in the literature. 

2. This turnover in the number counts is associated with 
a reddening of the faint ERO population, from a median 
colour of {hiA-K) = 4.38±0.05 at if ^ 19.0 to {hu-K) = 
4.71 ± 0.08 at 19.0^ K s;20.0. 

3. Our visual morphological fraction s at the bright end 
(K^19.0) are in good agreement with lYan fc ThomDsonI 
(■200S) . and we extend this sample to fainter depths, find- 
ing a mix of 35% disk-like, 15% disturbed/irregular, 30% 
spheroidal or compact and 20% unclassifiable. 

4. Using quantitative measures, we find that the concen- 
trations of EROs evolves across the break in their number 
counts, such that the fraction of disk-like galaxies declines. 

5. W e find on the basis of I / J / K colours, using the 
scheme of IPozzetti fc Mannuccil i2000f) . that most disk-like 
and amorphous galaxies are associated with dusty star- 
formation activity, and most bulge-like and compact sources 
have the colours of passive stellar populations. However, 
~30% of disks appear to have the colours of passive stel- 
lar populations. These are probably genuinely passive disks 
or at least have substantial evolved stellar populations. 
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6. We also fit the SEDs of our EROs finding sligfitly bet- 
ter agreement between tlie morpliological classification and 
the expected st ar-formation activity of such morphological 
types, than the IPozzetti fc Mannucci ll200(f) classification, 
i.e. most disk-like EROs exhibit photometric signs of recent 
star-formation. 

7. For a subsample of our EROs we derive SED-fitted 
photometric redshift distributions with a median of z = 
1.20 ±0.22, in reasonable agreement with spectroscopic red- 
shifts of small samples of somewhat brighter EROs. 

We have highlighted the diversity of Extremely Red Ob- 
jects, comparing morphological and photometric classifica- 
tion schemes. We suggest that the apparent break in ERO 
counts at K~19 is due to a rapid rise in bluer disk EROs 
at bright magnitudes. The dominance of this population de- 
clines as it is joined around K ^ 18.5 by a more modestly 
increasing population of galaxies comprising bulges, amor- 
phous, merging and increasingly unclassifiable systems. We 
have shown that these blue disk EROs have properties not 
inconsistent with them being at the same redshifts as the 
rest of the other ERO sub-populations. 

The next steps in understanding the nature of the ERO 
population are: spectroscopic follow up of a large samp le of 
faint EROs (building on the initial work of iCimatti et alJ 
|2002| at brighter magnitudes) in order to obtain star for- 
mation rates, extinction values and redshift distributions 
(particularly to understand the bluer disk-like EROs); and 
a detailed comparison with semi-analytic and other galaxy 
formation models. 
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Figure 11. Thumbnail images of a examples of EROs. Each set of images comprises an F814W and ii's-band pair; ERO ID number is 
given in upper left. Left panel (labelled with R.A. and Dec.) is the non-resampled F814W image, rebinned 2x2 in its original orientation; 
right panel is Xg-band image, rotated to match the F814W orientation, with the vector indicating the orientation on sky. The label gives 
the (7 — if) colour, K magnitude and visual morphological class in square brackets. The contours show the i^s-band isophotes. EROll is 
the largest disk in our sample. The Xs-band images appear overly smoothed as they have been resampled to match the (non-resamplcd) 
i814 image in this illustration, but are not resampled in our analysis. 



Table 6. Example of our ERO catalogue. Columns show: ERO ID; Field name (derived from MDS naming); RA; Dec; Total if-band 
magnitude; I-K colour, MDS Maximum Likelihood classification; concentration index; visual morphological class; and comments. The 
full table is available in the electronic version of the journal. 



ID# Field a (J2000) 5 (J2000) K-Tot {I - K) MLE C Vis Comments 



001 


U2AY2 


15:58:48.7 


+42:05:14 


19.69+0.14 


4.60+ 0.13 




0.30 


2 




002 


U2AY2 


15:58:50.5 


+42:05:03 
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4.33+ 0.13 


bulge 


0.25 
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early-type? 
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disk 
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or merger (near chip edge) 
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